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ABSTRACT 

The Red MSX Source (RMS) survey has identified a sample of ~1200 massive young stellar 
objects (MYSOs), compact and ultra compact Hii regions from a sample of ~2000 MSX and 
2MASS colour selected sources. We have used the 100-m Green Bank telescope to search 
for 22-24GHz water maser and ammonia (1,1), (2,2) and (3,3) emission towards ~600 RMS 
sources located within the northern Galactic plane. We have identified 308 H2O masers which 
corresponds to an overall detection rate of ~50 per cent. We find no significant difference in the 
detection rate for Hii regions and MYSOs which would suggest that the conditions required 
to produce maser emission are equally likely in both phases. Comparing the detection rates 
as a function of luminosity we find the H2O detection rate has a positive dependence on the 
source luminosity, with the detection rate increasing with increasing luminosity. 

We detect ammonia emission towards 479 of these massive young stars, which corre- 
sponds to ~80 per cent. Ammonia is an excellent probe of high density gas allowing us to 
measure key parameters such as gas temperatures, opacities, and column densities, as well as 
providing an insight into the gas kinematics. The average kinetic temperature, FWHM line 
width and total NH3 column density for the sample are approximately 22 K, 2kms ' and 
2 X 10'^ cm^"^, respectively. We find that the NH3 (1,1) line width and kinetic temperature 
are correlated with luminosity and finding no underlying dependence of these parameters on 
the evolutionary phase of the embedded sources, we conclude that the observed trends in the 
derived parameters are more likely to be due to the energy output of the central source and/or 
the line width-clump mass relationship. 

The velocities of the peak H2O masers and the NH3 emission are in excellent agreement 
with each other, which would strongly suggest an association between the dense gas and the 
maser emission. Moreover, we find the bolometric luminosity of the embedded source and the 
isotropic luminosity of the H2O maser are also correlated. We conclude from the correlations 
of the cloud and water maser velocities and the bolometric and maser luminosity that there 
is a strong dynamical relationship between the embedded young massive star and the H2O 
maser 

Key words: Stars: formation - Stars: early-type - ISM: molecules - ISM: radio lines. 



1 INTRODUCTION 



* -r ur T J c J f rr • f c- :i j -7 1 1 ui Masslve vounE Stellar objects (MYSOs) and ultra compact (UC) 

tables 2 and 5 and mil versions 01 rigs. 3 and 7 are only available j c jv / f \ / 

in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ^11 regions are two of the earliest phases in the lives of OB stars. 

(130.79.125.5) or via http://cds vyeb.u- st rasbg.fr/cgi-bin/qcat?J/A+A/ The MYSO phase begins when heating of the proto stellar enve- 

t E-mail: James.Urquhart@csiro.au (CASS) lope increases visibility of the object in the mid-infrared, and ends 
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once the central star begins to ionize its surrounding environment 
and forms an UC Hn region. These two phases are physically dis- 
tinct from the earlier hot molecular core phase, wh ich is not gener- 
ally d etectable at mid-infrared wavelengths (e.g.. lOe Buizer et alj 
2002h. MYSOs also possess strong ionized stellar winds (e.g., 



Table 1. Observed transition frequencies and excitation temperatures. 



Bunn et al.l 19951) : however, the radio emission from these winds is 



relatively weak (~1 mjy at 1 kpc uHoare et al.l 19941) and easily dis 
tinguishable from the slightly later radio-loud UC Hn region phase. 
It is likely that many MYSOs and UC Hn regions are still accreting, 
as evidenced by their almost ubiquitous association with powerful 
bipolar outflows JLadall 19851) . 

Models suggest that high levels of accretion in the early stages 
of a massive star's development lead to a 's welling up' of the pro- 
tosteUar core due to trapped entropy (e.g., Yorke & Bodenheimen 
l2008l : iHosokawa & Omukail |2009| ; iHosokawa et all I2OI0I) . The 
length of this swollen phase is comparable to the Kelvin-Helmholz 
timescale, which for stars with a current mass of >20 Mq is only 
a few thousand years. This means that once the accreted mass ex- 
ceeds ~20 Mq the star arrives on the main-sequence , even if accre- 
tion is still ongoing (see also lMcKee & Tanll2003h . At this point, 
the star begins to ionize its surroundings and becomes a UC Hn re- 
gion. This means that the MYSO phase is likely to be very brief 
(~ 10^ yr s) and is limited to objects with current masses lower than 
~20Mo l IDavies et ahluOUh . The UCHii region phase lasts a fac- 
tor of 2-5 times longer, with a small dependen ce on Ihe final stellar 
mass JDavies et al.ll201 itJMottram et al]|201 ih. 

The Red MSX Source (RMS; iHoare et alj |2005| ; 
iMottram et"ai] l2006l : lUrguhart era!] |2008|) Survey has estab- 
lished a large (~1200) and well-selected sample of MYSOs and 
compact and UCHn regions, and a database of complementary 
multi-wavelength datajj We have used arcsecond-resolution 
mid-infrared imaging from the Spitzer GLIMPSE survey 



I Benjamin & et al. 20031) or our own ground-based imaging (e.g., 
iMottram et alJuOOTH to reveal multiple and/or extended sources 
within the MSX beam, as well as MYSOs in close proximity to 
existing Hii regions. We have obtaine d arcsecond-resolution radio 
continuum with ATCA and the VLA JUrquhart et"al]|200l[200^ 
to identify UC Hn regions and PNe, w hilst observations of '^CO 
transitions JUrquhart et alj2007ll2008h deliver kinematic distances 
and luminosities, which allow us to distinguish between nearby 
low- and intermediate-mass YSOs and genuine MY SOs. Finally 
we h ave obtained near-infrared spectroscopy (e.g., Iciarke et al.l 
120061) which allows us to distinguish the more evolved stars. 

With the source classification effectively complete, the next 
step is to examine the global characteristics of this Galaxy-wide 
sample of massive young stars. This involves determining the phys- 
ical and chemical nature of the environment as a way of gauging the 
evolutionary status of our sample of MYSOs and Hn regions. In this 
paper we present the results of a set of thermal ammonia (NH3) and 
water (H2O) maser observations made with the 100-m Green Bank 
telescope (GBT) towards a sample of ~600 young massive stars. 

Ammo nia is an excellent tracer of high-density gas 
(-10** cm-^: lEvanslll999l : IStahler & Palldl2005l) . and is relatively 
unaffected by depletion at lower temper atures compared to othe r 
common molecular tracers such as CO JBergin & Langenll997n . 
The NH3 (1,1), (2,2) and (3,3) inversion transitions are normally 
coUisionally excited and the rotational temperature of the gas can 
be determined from the intensity ratio of any two inversion tran- 
sitions. The inversion transition splits into 18 separate hyperfine 



www.ast.leeds.ac.uk/cgi-bin/RMS/RMS.DATABASE.cgi. 



Transition 



Frequency" Energy Maser"^ Sensitivity 

(GHz) (K) or Thermal 



NH3(1,1) 23.694471 

NH3 (2,2) 23.722634 

NH3 (3,3) 23.870130 

H2O (61.6 -52,3) 22.235180 



23.4 


thermal 


52 mK 


64.9 


thermal 


50 mK 


124.5 


both 


54 mK 




maser 


0.12Jy 



" Transition rest frequencies have been taken from the Lovas catalogue 

available from|http://physics. nist.gov/cgi-bin/micro/table5/start.pI I 

'' Ener gy above the ground state in equivalent temperatures taken from 

lHo&T ownesll98J). 

^ Indicates whether a transition is a thermal line, a maser or possibly both. 



components JHo & Townejl 19831) , usually resolved into five dis- 
tinct components, the ratio of which can be used to calculate the 
optical depth of the transition. 

H2O masers are kno wn to occur in both high and low-mass 
star-f orming regions (e.g.. lForster & Caswelllll999l : IClaussen et al.l 
1 1 9961) and are thus an important signpost of ongoing star formation. 
Water mas ers are gener ally thought to be associated with molecular 
outflows ( iCodella et al.li2004l and references therein). We will use 
the H2O detection rates for our sample of MYSOs and Hn regions 
to investigate their statistical association as a function of evolution- 
ary phase and luminosity. 

The structure of the paper is as follows: in Sect. 2 we discuss 
the source selection, the observational setup and the data reduc- 
tion procedure. In Sect. 3 we use the reduced spectra to determine 
physical properties of the star-forming environments. In Sect. 4 we 
present the results and discuss their implications. We present a sum- 
mary of our results and highlight our main findings in Sect. 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Source selection 

The multi-wavelength data sets compiled as part of the RMS sur- 
vey have been used to cla ssify the initial sampl e into a number of 
different source types (see lUrquhart et alj|2008l for an overview of 
the classification scheme). For these observations we have selected 
all RMS sources classified as a young stellar object or Hn region. In 
some cases there is more than one embedded source located within 
the MSX beam (~18") or a source can display traits of both Hn 
regions and YSOs and is possibly in a transitional stage between 
the two: these sources are classified as 'Hn/YSO' and were also in- 
cluded in our sample. Finally, we have been unable to definitively 
classify approximately 10 per cent of sources where the available 
data are contradictory; these are designated 'young/old' to indi- 
cate the uncertainty associated with them. Although the majority 
of these sources are likely to be evolved stars, we include them to 
avoid exclusion of any genuine YSOs that may still be lurking in 
this group. 

This selection process produced a sample of 586 sources lo- 
cated in the northern Galactic plane (i.e., 10° < 1 < 180°). This 
sample provides a complete census of compact and UC Hn regions, 
and MYSOs in the first and second quadrants. We complement this 
representative sample with a further eleven sources located in the 
outer part of the southern Galactic plane (i.e., 220° < I < 240°). 
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Figure 1. Plots of the noise distribution of tlie four observed lines. 



2.2 GBT Observations 

Observations were made of the NH3 (1,1), (2,2) and (3,3) inversion 
transitions and the H2O maser transition towards a sample of 597 
young massive stars identified by the RMS survey (see Table[T]for 
transition rest frequencies). These observations were made during 
six sessions starting on the 25* of November 2009 and ending on 
the 10* of December 2010 using the Green Bank telescope (GBT), 
operated by the National Radio Astronomy ObservatorjtJ. 

All three NH3 transitions were observed simultaneously, elim- 
inating many of the sources of uncertainty while allowing many 
key parameters such as excitation, rotation and kinetic tempera- 
tures, column densities, and optical depths to be calculated for a 
large sample of MYSOs and Hii regions. The observations were 
performed in frequency-switch mode in order to remove sky contri- 
butions and a noise diode was observed in switching mode through- 
out the observations in order to achieve absolute flux calibration on 
the r^ scale to an accuracy of ~ 10 per cent. The natural spectral 
resolution of the raw data was 6. 1 kHz (~0.08 km s"'). A standard 
integration time of 70 seconds was used for each source, resulting 
in a typical r.m.s. noise of 100 mK per channel with median system 
temperatures of 61 K. Weather conditions were stable for all ob- 
serving sessions: typical pointing offsets were ~5" in azimuth and 
~3" in elevation with a half-power beam- width of ~30". 
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Figure 2. Ammonia spectra towards the MSX source G133. 9476+01.0648 
as an example of a source with both emission and absorption. The full ver- 
sion of this figure can be found in Fig lAll 



The reduced spectra were Hanning-smoothed to a resolution 
of ~0.32 km s"' and a sensitivity of 50 mK channel"' . In Fig.[T]we 
present histograms of the noise distributions of the four observed 
transitions. Noise values are given in mK for the ammonia transi- 
tions; however, for the water masers we have converted the temper- 
ature scale to the more commonly used Jansky scale using the GBT 
gain value of 1.5 K Jy '0 

The final reduced spectra were calibrated to the corrected 
antenna temperature scale (T^), and converted to the telescope- 
independent main-beam temperature scale (r„,b), assuming a main- 
beam efficiency (i/mb) of 0.89'*. 

The NH3 emission seen towards the vast majority of sources 
can be attributed to a single molecular cloud located along the line 
of sight. Multiple clouds, identified through multiple spectral com- 
ponents, are seen towards only nine sources. More complex emis- 
sion structure is seen towards six where the spectra appear to be 
combinations of emission and absorption; these spectra indicate 
a bright continuum source in the beam (i.e.. Hit region). In all of 
these cases the velocity ranges of the absorption and emission over- 
lap, making it difficult to obtain reliable parameters. We present the 
emission detected towards one of these sources in Fig.[2](the spec- 
tra seen towards all six of these sources can be found in Fig. lAll of 
the Appendix). Since the results obtained by fitting these spectra 
are unlikely to be reliable, we simply present these spectra without 
further analysis. 



2.3 Data reduction 

The data were reduced using the GBTIDLJj data analysis package. 
Bad scans were removed and channels outside the region of inter- 
est were discarded. Temperature scale corrections for atmospheric 
opacity were made using the zenith values provided from local 
weather models. A high-order polynomial was fitted to emission 
free channels and subtracted from all spectra to remove baseline 
anomalies before any emission features were fitted. 



3 DERIVING PHYSICAL PROPERTIES 

In this section we will describe the methods used to determine 
physical properties of the dense environments in which young mas- 
sive stars are forming. 



^ The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
' http://gbtidl.nrao.edu/ 



"* The gain factor and main beam efficiency have been taken from 'The 
Proposer's Guide for the Green Bank Telescope'. 



Table 2. Observed parameters. 



NH3 (1,1) 



NH3 (2,2) 



NH3 (3,3) 



MSX Name" 


Type'' 




RA 


Dec 


Ta 


^roi 


T-lun 


Log[(V(NH3)] 


Tmb' 


Vlsr 


AV 


Tlni.l.l) 


Tmb 


Vlsr 


AV 


Tmb 


Vlsr 


AV 


H2O 






(J2000) 


(J2000) 


(K) 


(K) 


(K) 


(cm"-) 


(K) 


(kras"') 


(kms"') 




(K) 


(kms"') 


(kms"' 


) (K) 


(kms"') 


(kms"') 




GOIO.3844+02.2128 


1 


18 


00:22.68 


-18:52:08.0 


4.5 


14.8 


16.0 


15.0 


1.2 


5.54 


0.79 


1.11 


04 


5.44 


0.77 








y 


GOlO.4413+00.0101 


2 


18 


08:38.23 


-19:53:574 


3.7 


16.8 


18.5 


15.7 


0.6 


66.37 


2.63 


1.54 


0.3 


66.88 


3.08 


0.2 


66.28 


4.54 


n 


GOlO.4718+00.0206 


2 


18 


08:39.67 


-19:52:03.0 


4.1 


20.7 


24.0 


16.2 


1.0 


67.40 


5.73 


2.17 


0.8 


67.55 


5.90 


0.8 


67.50 


7.88 


y 


GOlO.5067+02.2285 


1 


18 


00:34.58 


-18:45:17.6 


5.7 


14.5 


15.5 


15.4 


24 


21.79 


1.31 


1.83 


0.9 


21.86 


1.69 


0.2 


21.94 


3.24 


y 


GOlO.6291-00.3385 


2 


18 


10:19.32 


-19:54:12.9 


6.9 


17.2 


19.1 


15.9 


3.6 


-4.54 


1.86 


3.39 


2.5 


-4.51 


2.57 


1.1 


-4.42 


3.19 


n 


G010.631 1-00.3864 


2 


18 


10:30.26 


-19:55:30.0 


15.1 


27.8 


35.5 


15.4 


2.0 


-3.49 


5.07 


0.25 


1.6 


-4.84 


6.41 


2.0 


-4.95 


7.56 


y 


G010.841 1-02.5919 


1 


18 


19:12.16 


-20:47:32.2 


6.1 


23.7 


28.5 


15.1 


1.0 


12.27 


1.91 


0.41 


0.7 


12.25 


2.16 


04 


12.33 


2.36 


y 


GOlO.8856+00.1221 


1 


18 


09:08.13 


-19:27:23.0 


4.0 


20.1 


23.0 


15.3 


0.6 


19.67 


2.13 


0.77 


0.3 


19.69 


2.51 


0.1 


19.84 


3.70 


y 


GOlO.9592+00.0217 


2 


18 


09:39.57 


-19:26:25.8 


6.5 


93.8 


28.8 


15.1 


0.9 


2141 


2.90 


0.29 


0.6 


21.61 


3.51 


04 


21.20 


4.24 


y 


GOlO.9657+00.0083* 


2 


18 


09:43.36 


-19:26:28.6 










0.2 


19.16 


3.98 




0.2 


18.36 


3.20 


0.1 


19.42 


1.87 


y 


GOll. 1109-00.4001 


2 


18 


11:32.27 


-19:30:40.6 


5.1 


19.1 


21.6 


15.7 


1.6 


0.12 


3.04 


1.33 


0.9 


0.01 


3.68 


0.5 


-0.14 


4.30 


n 


GOU. 1723-00.0656 


2 


18 


10:25.22 


-19:17:46.3 










<0.2 




















n 


GOl 1.3252-01.8040* 


4 


18 


17:12.89 


-19:59:36.9 










04 


9.00 


0.95 
















y 


GOl 1.3757-01.6770* 


2 


18 


16:50.47 


-19:53:204 










0.2 


14.62 


2.06 




0.1 


14.49 


1.74 








n 


GOl 1.3757-01.6770* 


2 


18 


16:50.47 


-19:53:204 










0.3 


7.68 


1.01 




0.2 


7.76 


0.96 








n 


GOll.4201-01.6815 


3 


18 


16:56.86 


-19:51:07.2 


6.0 


17.8 


19.9 


15.5 


2.5 


9.18 


1.88 


1.50 


1.3 


9.20 


2.22 


04 


9.23 


2.56 


n 


GOll.5001-014857 


3 


18 


16:22.58 


-19:41:19.3 


8.2 


214 


24.9 


14.9 


14 


10.40 


1.74 


0.33 


0.8 


1041 


1.97 


0.2 


10.24 


3.26 


y 


GOl 1.9019+00.7265* 


1 


18 


08:58.89 


-18:16:28.2 










0.2 


24.28 


1.77 




0.2 


24.46 


1.69 








n 


GOl 1.9454-00.0373 


2 


18 


11:53.47 


-18:36:18.3 










<0.2 




















n 


G012.0260-00.0317 


1 


18 


12:02.04 


-18:31:54.1 


3.7 


18.7 


21.1 


15.6 


0.7 


110.85 


1.97 


1.59 


04 


110.81 


2.50 


0.2 


110.56 


341 


n 


G012. 1993-00.0342 


3 


18 


12:23.66 


-18:22:51.6 


4.6 


19.5 


22.2 


15.7 


1.3 


50.86 


2.62 


1.52 


0.8 


50.88 


3.21 


0.5 


51.65 


4.71 


y 


G012.4314-01.1117 


2 


18 


16:51.33 


-18:41:30.1 


5.9 


214 


24.9 


15.4 


1.9 


39.73 


1.97 


0.99 


1.2 


39.68 


2.44 


0.9 


39.57 


3.20 


y 


G012.5932-00.5708* 


4 


18 


15:10.43 


-18:17:30.8 










0.2 


15.87 


1.60 
















n 


G012.8600-00.2737 


2 


18 


14:36.67 


-17:54:56.5 


5.7 


17.8 


19.9 


15.8 


2.3 


36.77 


2.15 


2.34 


14 


36.69 


2.65 


0.6 


36.60 


2.76 


n 


G012.8909+00.4938 


3 


18 


11:50.69 


-17:31:14.5 


6.1 


20.1 


23.0 


16.0 


2.3 


32.92 


3.22 


2.23 


1.6 


32.83 


3.87 


1.0 


32.81 


4.24 


y 


G012.9090-00.2607 


1 


18 


14:39.71 


-17:51:59.0 


8.5 


21.9 


25.7 


15.9 


4.0 


36.73 


3.54 


1.55 


2.8 


36.88 


3.68 


2.0 


36.89 


4.55 


y 


G013.0105-00.1267 


4 


18 


14:22.25 


-17:42:47.8 


3.8 








0.4 


10.42 


3.23 


0.85 














n 


G013.0105-00.1267 


4 


18 


14-99 25 


-17:42:47.8 


3.8 








04 


10.43 


3.20 


0.82 














y 


GOB. 1840-00.1069 


3 


18 


14:38.76 


-17:33:05.0 


5.2 


16.6 


18.2 


15.7 


2.1 


53.39 


1.61 


2.83 


1.2 


53.38 


2.07 


0.6 


53.18 


3.60 


n 


GOB. 1840-00.1069 


3 


18 


14:38.76 


-17:33:05.0 


4.1 








1.0 


36.79 


0.86 


1.45 














n 


G013.3310-00.0407 


1 


18 


14:41.78 


-17:23:26.1 


5.9 


154 


16.7 


15.6 


2.7 


54.68 


1.44 


2.39 


1.3 


54.65 


1.87 


04 


54.51 


3.79 


n 


G013.6562-00.5997 


1 


18 


17:24.36 


-17:22:14.5 


6.0 


18.2 


20.4 


15.9 


2.8 


47.45 


2.08 


3.52 


2.1 


47.49 


2.84 


2.0 


47.58 


3.31 


y 


G013.8885-00.4760 


3 


18 


17:24.67 


-17:06:27.3 


5.6 








0.7 


22.68 


0.92 


0.31 














n 


G014.0329-00.5I55 


1 


18 


17:50.61 


-16:59:57.1 


5.7 


20.9 


24.3 


14.8 


0.6 


20.34 


1.47 


0.31 


0.3 


20.42 


1.39 


0.1 


20.63 


2.12 


n 


G014.2166-00.6344 


1 


18 


18:38.70 


-16:53:37.3 


3.8 


14.6 


15.7 


15.2 


0.6 


19.75 


1.36 


1.10 


0.2 


19.74 


1.78 








y 


G014.3313-00.6397 


2 


18 


18:53.49 


-16:47:42.7 


7.2 


20.8 


24.1 


15.8 


2.8 


22.26 


9.87 


1.55 


1.9 


22.31 


3.23 


1.5 


22.65 


3.55 


y 


G014.4335-00.6969 


1 


18 


19:18.23 


-16:43:55.9 


5.2 


18.6 


21.0 


15.1 


1.2 


17.50 


1.38 


0.76 


0.6 


17.65 


1.66 


0.2 


17.96 


1.97 


n 


G0144886+00.0219 


2 


18 


16:46.27 


-16:20:344 


3.5 


17.3 


19.1 


15.8 


0.5 


22.88 


3.01 


1.70 


0.3 


22.90 


4.05 


0.1 


22.69 


4.13 


n 


G014.5982+00.0202 


2 


18 


16:59.68 


-16:14:50.6 


4.1 


19.9 


22.8 


15.6 


1.0 


26.62 


1.90 


1.48 


0.6 


26.56 


2.53 


04 


26.72 


4.64 


y 


G014.6087+00.0I27 


1 


18 


17:02.56 


-16:14:30.1 


5.1 


20.1 


23.1 


15.9 


1.8 


24.63 


2.92 


2.02 


1.3 


24.63 


3.64 


1.1 


24.60 


4.50 


y 


G014.9790-00.6649 


3 


18 


20:15.60 


-16:14:10.3 


3.7 


20.9 


24.2 


15.0 


0.6 


18.91 


0.76 


1.04 


0.3 


18.83 


1.17 


0.2 


19.09 


1.70 


n 


G014.9958-00.6732 


1 


18 


20:19.43 


-16:13:31.0 


4.8 


23.1 


27.7 


16.0 


1.1 


19.43 


4.79 


1.34 


0.8 


19.22 


5.28 


0.5 


19.07 


5.04 


y 


G015.0357-00.6795* 


2 


18 


20:25.51 


-16:11:35.5 










04 


19.20 


2.59 




0.5 


19.40 


2.85 


0.6 


19.70 


2.78 


y 


G015.0939+00.19B* 


1 


18 


17:20.86 


-15:43:47.2 










0.1 


29.93 


2.02 
















y 


G015. 1288-00.67 17 


1 


18 


20:34.75 


-16:06:26.2 


3.1 


24.8 


30.3 


15.5 


0.3 


18.99 


1.17 


1.63 


0.2 


18.80 


1.74 


0.1 


17.67 


343 


y 


G016.1438+00.0074 


2 


18 


20:04.77 


-14:53:31.2 


3.9 


16.2 


17.7 


15.5 


0.9 


44.84 


1.84 


1.54 


04 


44.86 


2.46 


0.2 


44.90 


3.75 


y 


G016.7122+01.3119 


1 


18 


16:26.85 


-13:46:24.9 


5.8 


15.2 


16.5 


15.1 


1.8 


20.31 


1.12 


0.93 


0.6 


20.26 


1.39 








n 


G016.798I+00.1264 


1 


18 


20:55.27 


-14:15:32.0 


4.2 


14.1 


15.1 


15.3 


1.2 


15.19 


0.71 


2.44 


0.5 


15.19 


0.88 








n 


G016.8055+00.8149 


1 


18 


18:25.96 


-13:55:37.9 


5.5 


13.0 


13.7 


14.9 


2.0 


19.73 


0.54 


1.38 


0.5 


19.75 


0.70 








y 


G016.8689-02.1552 


1 


18 


29:24.26 


-15:15:42.1 


8.0 


18.8 


21.2 


16.0 


44 


18.58 


2.99 


2.44 


3.0 


18.54 


3.82 


2.0 


18.49 


5.03 


y 


" A superscript dagger symbol appended to the MSX name indicate source detections where we have been uirable to fit the hyperfine satellite features of the NH3 (1,1) emission 
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NH^ and H2O maser analysis of massive star forming regions 
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Figure 3. Example of the NH3 emission detected towards a young mas- 
sive star, with NH3 (1,1), (2,2) and (3,3) spectra on the top, lower 
left and lower right plots respectively. The model fits to these data are 
shown in red. The full version of this figure includes spectra taken 
towards all 591 sources and is only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-i-A/] 



3.1 NH3 : Molecular gas properties 

An example of the reduced spectra of the NH3 (1,1), (2,2) and (3,3) 
inversion transitions and the fit to the data is presented in Fig.[3] 
We detect NH3 (1,1) emission towards ~80 of the sample with the 
hyperfine structure seen towards the majority (~65 percent) of the 
observed sources, and used an IDL routine to simultaneously fit 
all 18 hyperfine components and to derive the optical depth and 
line widths. The hyperfine structure is generally too weak to be ob- 
served in the NH3 (2,2) and (3,3) transitions: for these lines (as well 
as the remaining NH3 (1,1) lines) we have obtained corrected an- 
tenna temperatures by fitting a single Gaussian profile to the main 
line. In all cases the NH3 (1,1) and (2,2) line widths have been ob- 
tained by fitting the hyperfine components to their respective main 
line emission to remove the effects of line broadening due to opti- 
cal depth. The resulting fits to the data are shown in Fig.[3]and the 
fitted parameters, velocity, line width, optical depth and corrected 
antenna temperature are presented in Table l2l^l 

We obtain the total optical depth (T(i_i)) of the NH3 (1,1) tran- 
sition from the hyperfine fitting performed. If we assume that the 
optical depth ratios between the hyperfine transitio ns are equal to 
the ra tios of the theoretical transition intensities JRvdbeck et al.l 
ll977D , then the optical depth of the NH3 (1,1) main quadmpole 
transition (T(i !_,„)) is approximately half that of the total optical 
depth. We can now estimate the excitation temperature {T^^) using 
the detection equation assuming the gas is in local thermal equilib- 
rium (LTE), a background temperature of 2.73 K and a beam filling 
factor (^ff ) of unity; the beam filling factor is a measure of the frac- 
tion of the telescope beam that is filled by the source. Using the 
main line optical depth (T(i. !_,„)) of the NH3 (1,1) transition, we ob- 
tain the excitation temperature using: 
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Figure 4. Distributions of measured and derived parameters for the ob- 
served sample. The mean and standard derivation, and maximum and mini- 
mum values of each distribution are presented in Table|3] 



n = '7mb'7ff[^(re.) - /(7'bg)][l - e-^"'-']. 



where. 



UT) = 



hv/k 



(eh./kr _ 1) • 



(1) 



(2) 



The rotation temperature (7^1) associated with the NH3 (2,2) 
and (1,1) transitions can be calculated using the line intensities and 
the (1 ,1) main quadrupole transition optical depth JHo & Townej 
Il983h : 



^ The full version of Fig.|3]and Table|2]aie only available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-l-A/| 



In 



-0.282 
r, l.l.m) ' 



In l-# 



■d 



_e-'"(i.i,m)) I 



[K], (3) 



where To = '^~'^, "'" » 41.5 K is the temperature associated with 
the energy difference between (1,1) and (2,2) levels, and ^^^ , ^^ 
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Table 3. Summary of derived parameters 



Pai'ameter 


Number 


Mean 


Median 


lo- 


Min 


Max 


Tmb (K) 


479 


1.1 


0.7 


1.0 


0.1 


6.9 


AV(kms-l) 


479 


1.9 


1.7 


0.9 


0.2 


7.8 


m.i) 


386 


2.2 


2.1 


1.6 


0.3 


9.2 


T-ex (K) 


386 


5.3 


4.9 


1.5 


3.1 


15.1 


Trot (K) 


365 


19.3 


19.2 


3.8 


10.2 


39.0 


Tkm (K) 


365 


22.1 


21.4 


5.9 


10.5 


60.8 


Log[N(NH3)(cm-2)] 


365 


15.3 


15.3 


0.4 


14.4 


16.5 


Beamjf 


365 


0.28 


0.27 


0.09 


0.09 


0.63 


Log[Sy(H,0) (Jy)] 


298 


1.1 


0.9 


1.0 


-0.5 


5.8 


Log[j5,,(H20) (Jykms"')] 


298 


1.3 


1.2 


1.1 


-1.7 


6.7 


H20Vel. Range (km s-') 


298 


24.7 


15.2 


32.5 


0.00 


362.0 


Log[L(H20) (Lo)] 


298 


-5.0 


-5.0 


1.3 


-8.9 


1.2 



and 7^(^2111) ^^ ^^ measured corrected peak antenna tempera- 
tures of the main quadrupole transition of the NH3 (1,1) and (2,2) 
hues, respectively. The next step is to relate the rotation tempera- 
ture to the kinetic temperature of the gas. Empirical results reveal 
that at low temperatures (< 15 K) the rotation and kinetic temper- 
ature are approximately equivalent; however, they begin to devi- 
ate at higher temperature and thus analytic expressions underesti- 
mate the rotation temperatures for kinetic temperatur es above 40 K 
jWalmslev & Ungerechts|[T983l : lHo & Towneslll983l) . For T,, < To, 
a relationship between excitation and kinetic temperature may be 
calculated by consideration of (1,1), (2,2), and (2,1) states only 
JSwift et al.ll2005l : IWalmslev & Ungereditsll 19831) . such that: 



^rot — 



Tkin 



1 + 



^ln[l+0.6xexp(^)] 



[K] 



(4) 



The kinetic temperatures calculated using Eqn. 3 may be 
overestimated for sources where r^ot - Tq. The mean rotation 
temperature is ~20K and so this is not a concern for the vast 
majority of our detections. However, we do find a handful of 
sources where this might be an issue (e.g., G019.6085-00.2357 
and G049.4903-00.3694). 

If excitation conditions are homogeneous along the beam and 
all hyperfine lines have the same excitation temperature, then the 
column densities at a given {J, K = J ) transition can be wr itten as a 
function of the total column density JMangum et al.ll 19920 : 



N, 



(1,1) 



6.60 X 10"AV(,.i)T( 



i,i,m) [cm ] 

m.i) 



(5) 



where A'(i.i) is the column density of the NH3 (1,1) transition, 
AV(i,i) is the FWHM line width of the NH3 (1,1) transition in 
km s"' £ind V(i,i) is the transition frequency in GHz . Finally we est i- 
mate the total ammonia column density following lLi et alj ( 120031) : 



M 



NH, 



^A'd 



, 1 /23.1 



+ -exp 



-41.2 



14 
+ -exp 



-99.4 



(6) 
In Fig.|4]we present histograms showing the distributions of 
the various parameters derived in this subsection and present a sum- 
mary of the derived parameters in Table[3] 
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Figure S. Plot comparing the derived excitation and rotation temperatures. 
The red dashed line indicates the line of equality. The black vertical line in- 
dicates the minimum cloud temperatures expected due to heating by cosmic 
rays and the far-uv radiation field (~10 K: lEvansl 19991) . whilst the black hor- 
izontal line indicates the temperature of the cosmic microwave background 
radiation. 



3.2 Beam filling factor 

In Eqn. 1 presented in the previous subsection we assumed the 
beam filling factor to be of order unity; however, having calculated 
the excitation and rotation temperatures of the molecular gas we 
are in a position to test this assumption. In Fig.[5]we present a scat- 
ter plot comparing the derived excitation and rotation temperatures; 
the dashed line indicates the line of equality where T^^ = Ty^j. This 
plot clearly shows that the rotation temperatures are systematically 
~4 times higher than the excitation temperatures. The excitation 
temperatures are calculated from a single transition and assume the 
beam is uniformly filled (i.e., rjg ~ 1). The assumption that the 
beam is uniformly filled results in unfeasibly low inversion transi- 
tion excitation temperatures (typically ~5 K). The rotation temper- 
ature is derived from the line intensity ratio of the NH3 (1,1) and 
(2,2) transitions and therefore the beam filling factor is effectively 
divided out and thus provides a more reliable estimate of the gas 
temperature. 

The assumption of LTE includes the inversion transitions and 
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Figure 6. Distribution of the differences between the velocity assigned from 
CO data and the velocity of the NH3 (1,1) inversion transition. In the 
upper left corner we give the statistical parameters of the distribution; these 
have been calculated for all sources where the absolute velocity difference 
is less than 5kms"'. Note the j:-axis has been truncated ±7kms"'; the 
details of the two sources that have velocity differences greater than this 
can be found in Table|4] 



implies tiiat Tex siiould be equal to r,ot. Hence, our derived values of 
Tex are actually ^ r,-o, x rjg and we can estimate the beam filling fac- 
tor from the ratio of the excitation and rotation temperatures, i.e., 
'7 If = T^^/Trof We note the non-linear nature of Eqn. 1, however, 
given that J{T} ~ T in Eqn. 2 for T » hv/k, which corresponds 
to a temperature of 1 . 1 K at v= 23 GHz and so is nearly always the 
case for these transitions, we find that Eqn. 1 is approximately lin- 
ear in T. We derive minimum and maximum beam filling factors 
of 0.1 and 0.6 respectively, with a mean value of ~ 0.3; this com- 
pare well to the values reported bv lPillai et alj ( l200q) (rj^ ~0.3-0.5) 
from ammonia mapping observations of a small sa mple of infrared 
dark clouds (IRDCs) and lRosolowskv et al.l ( I2OO8I) (77,, -0.3) from 
GBT observation s of dense cores in Perseus. The observations of 
IPillai et alj ( l200q) were conducted with the Effelsberg 100-m tele- 
scope and have a similar resolution and sensitivity to the obser- 
vations presented here. Since it is unlikely that the inversion transi- 
tions would be sub-thermally excited, whilst the rotation transitions 
are thermalised, the implication is that the low beam filling factors 
are real and so there is likely to be high-contrast substructure within 
the GBT beam. 

Further evidence of the presence of substructure comes 
from high-resolution observations of a sample of high-mass star 
forming clumps obser ved with the Compact Array reported by 
iLongmore et alj 1 120071) . These authors found typical clump sizes of 
~10-15" with multiple cores being detected within the 2' primary 
beam in approximately half of the fields observed. Using the beam 
filling factor we estim ate the typical angular diameter to be simi- 
lar to those reported bv lLongmore et al. 12007), which at the mean 
distance of our sample (~5 kpc), corresponds to physical diameters 
of 0.2-0.4 pc. These observations are therefore probing molecular 
clumps that are likely going to form single stars or small multiple 
stellar systems rather than whole clusters. 
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Table 4. RMS sources where the diiferences between the assigned '^CO ve- 
locity and NH3 (1,1) velocities are greater than 3o" (i.e., |AV| > 2.1 kms"'). 

MSXName Vlsr("CO) Vlsr(NH3) ^VJ 

(kms"') (kms"') (kms"') 

G031.4134+00.3092 
GOlO.4413+00.0101 
G035.1992-01.7424 
GO 16.95 12+00.7806 
G045.4658+00.0457 
G014.4335-00.6969 
G030.7206-00.0826 
G015.1288-00.6717 
G038.2577-00.0733 
G030. 1043-00.0738 
G017.1141-00.1130 



3.3 NH3-CO velocity comparison 

The RMS project has previo usly used '^CO rotational tra nsitions to 
determine source velocities JUrguhart et al.ll2007l uOOSf) and kine- 
matic distances, using the Galactic rotation curve of lBrand & Blitzl 
1 119930 . However, we detected multiple velocity components in 
'^^CO towards approximately 60 per cent of our sources. In the ma- 
jority of these cases it was straightforward to identify the correct 
velocity component. However, there were a small number in which 
the velocity assignment was less reliable. With the NH3 observa- 
tions in hand we are in a position to check the previous velocity 
assignments and make corrections where necessary. In Fig.[6]we 
present a histogram showing the difference between the assigned 
'^^CO velocity and the NH3 (1,1) velocity. A cursory inspection 
of this plot reveals the velocity assigned from the CO spectra is 
in excellent agreement with the NH3 velocities. Indeed there are 
only eleven sources where the velocities differ by more than 3cr 
(~2.1 kms"'); this corresponds to approximately 2 per cent of the 
sample. In Table|4]we give the names and velocities for these eleven 

sources. 

Inspection of the individual '""CO profiles JUrquhart et al.l 

uOOa) for all sources in Table|4]with velocity differences less than 
5 km s"' reveals the presence of two or more blended components, 
or a single emission feature with uncharacteristically broad line 
width of ~5-7kms"', which probably indicates that a number of 
velocity components have been blended into a single feature. The 
differences in the CO and NH3 velocities for these sources can eas- 
ily be understood in terms of blended components along the line of 
sight. The remaining three sources have velocity differences greater 
than 5 km s"' ; inspection of their CO spectra reveals the presence of 
several distinct emission features and, given that NH3 is a tracer of 
dense gas, it is likely that the wrong CO component was assigned. 
We have used the NH3 velocity to determine the kinematic prop- 
erties of these three sources and have updated the RMS database 
accordingly. 



3.4 Water maser properties 

In Fig.|7]we present the water-maser spectrum observed towards 
the MSX source GO43.1679+OO.OO950This particular source is an 



^ The full version of this figure, showing all 308 spectra, is available 
in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.125.5) or via |http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A7] 
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Table 5. Water maser properties. 



MSX Name 



Type" 



(Jy beam ) 



(kms-') 



(kms"') 



* max 



(kms"') 



Vnh, 
(kms-') 



Log[5„] 

(Jy) 



Log[j5„d\/] 
(Jykms"') 



Log[L(H,0)] 

(Lo) 



LogLLbd]" 

(Lo) 



G229.5711+00.1525 


2 


0.10 


40.4 


60.9 


46.7 


52.8 


1.80 


2.08 


-4.07 


4.11 


G236.8158+01.9821 


2 


0.11 


37.3 


65.2 


52.0 


52.8 


1.86 


2.43 


-3.76 


3.80 


G038.3543-00.9519 


2 


0.09 


11.2 


13.0 


12.3 


16.9 


0.57 


0.44 


-7.09 


1.96 


G232.0766-02.2767 


2 


0.09 


41.4 


48.5 


42.0 


42.1 


0.69 


0.74 


-5.64 


3.71 


G220.4587-00.6081 


2 


0.09 


23.6 


33.6 


30.5 


29.6 


1.02 


1.04 


-5.64 


3.18 


G224.6075-01.0063 


2 


0.08 


27.9 


29.7 


28.7 


16.7 


0.41 


0.38 


-6.86 


3.19 


G()79.8855+02.5517 


2 


0.10 


-2.0 


20.2 


4.2 


5.9 


0.99 


1.43 


-5.75 


3.55 


G079.8749+01.1821 


1 


0.10 


-7.4 


3.1 


-4.2 


-4.3 


1.27 


1.25 


-5.93 


3.83 


G082.0333+02.3249 


1 


0.09 


-4.8 


8.0 


-0.8 


0.3 


1.20 


1.55 


-5.63 


3.51 


G083.7071+03.2817 


2 


0.09 


-7.1 


8.4 


-6.3 


-3.6 


0.59 


0.86 


-6.32 


3.75 


G080.8624+00.3827 


2 


0.09 


-9.6 


4.2 


-7.7 


-1.9 


1.36 


1.64 


-5.54 


3.75 


G08 1.8652+00.7800 


2 


0.12 


-28.5 


46.5 


10.6 


9.4 


3.22 


3.81 


-3.37 


3.70 


G08 1.8789+00.7822 


1 


0.11 


-13.2 


20.9 


10.6 


8.1 


2.73 


3.30 


-3.88 


4.49 


G081.7131+00.5792 


2 


0.10 


-7.8 


15.2 


1.4 


-3.6 


1.34 


1.72 


-5.45 


3.83 


G08 1.7220+00.5699 


1 


0.10 


-44.8 


33.8 


-5.5 


-2.8 


3.06 


3.21 


-3.97 


2.80 


G08 1.7522 +00.5906 


2 


0.10 


-5.6 


27.4 


0.6 


-4.0 


0.84 


1.21 


-5.97 


3.39 


G081.7133+00.5589 


1 


0.10 


-6.6 


7.7 


-0.1 


-3.8 


0.31 


0.85 


-6.33 


3.42 


G08 1.7624+00.59 16 


2 


0.11 


-14.8 


32.9 


-9.1 


-4.4 


0.51 


0.96 


-6.22 


3.00 


G08 1.6632+00.4651 


2 


0.10 


8.1 


22.5 


8.4 


19.3 


0.63 


0.52 


-6.66 


2.54 


G084. 1940+01. 4388 


2 


0.10 


-0.2 


41.0 


38.7 


-1.8 


0.56 


0.89 


-6.29 


3.67 


G085.4102 +00.0032 


3 


0.10 


-39.6 


-12.9 


-32.5 


-35.8 


1.84 


2.24 


-3.74 


4.50 


G084.9505-00.6910 


2 


0.09 


-34.5 


-29.6 


-33.9 


-34.9 


0.20 


0.18 


-5.81 


4.28 


G095.053 1+03.9724 


2 


0.11 


-91.9 


-82.5 


-88.4 


-84.1 


1.63 


1.78 


-3.79 


4.28 


G093.1610+01.8687 


2 


0.11 


-70.0 


-70.0 


-70.3 


-63.1 


0.72 


0.36 


-5.43 


4.17 


G097.5268+03.1837 


3 


0.12 


-111.5 


-52.6 


-76.1 


-70.0 


2.42 


3.11 


-2.65 


4.69 


G094.2615-00.4116 


2 


0.10 


-76.8 


-35.6 


-36.4 


-46.6 


0.77 


1.04 


-5.02 


4.10 


G094.4637-00.8043 


2 


0.10 


-48.6 


-43.5 


-46.2 


-44.6 


1.03 


1.16 


-4.90 


4.50 


G096.4353+01.3233 


3 


0.11 


-70.0 


-62.5 


-69.9 


-68.9 


-0.03 


0.07 


-5.68 


4.27 


G094.6028-01.7966 


2 


0.12 


-60.2 


-39.8 


-56.0 


-43.9 


2.43 


2.88 


-3.64 


4.37 


G095.0026-01.5779 


2 


0.10 


-46.6 


-38.4 


-45.8 


-40.6 


1.41 


1.39 


-5.13 


3.48 



" Source types are as follows: (1) Hii regions; (2) YSOs; (3) Hii/YSOs; (4) Young/Ol d stars. 

' Bolometric luminosities have been calculated using the integrated fluxes derived bv lMottram et alj 1201 ih . 

Notes: Only a small portion of the data is provided here, the full table is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 

(130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bi n/qcat ?J/A+A/| 




Vlsr (km s"') 



Figure 7. Example of H2O maser emission detected towards the MSX 
source G043. 1679+00.0095. The vertical dashed lines indicate the min- 
imum and maximum velocity ranges over which maser emission is 
detected. The full version of this figure includes spectra taken to- 
wards 597 sources and is only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/] 



Hn region associated with the W49A star- forming complex, and the 
water maser detected towards this region is the brightest we have 



found with a peak flux value of ~ 6 x 10' Jy. Maser spectra of- 
ten consist of a number of distinct emission peaks, usually referred 
to as maser spots, spread over a range of velocities. The vertical 
dashed lines in Fig.|7]show the minimum and maximum velocities 
over which maser spots are found above a 3 cr detection limit to- 
wards this object. 

For all of the maser emission detected we measure the peak 
flux and velocity of the brightest maser spot, the integrated flux 
and the minimum and maximum velocities over which the maser 
spots are distributed. These values are tabulated in Table[5] In the 
upper panels of Fig.[8] we present the peak and integrated flux- 
density distributions and in the lower left panel of this figure we 
present a histogram of the total velocity range over which the maser 
spots are seen in each spectrum. The median peak flux is ~8Jy 
with minimum and maximum peak fluxes of 0.26 and 5.7x10' Jy 
respectively. The largest velocity range of maser spots is seen 
towards G043. 1679+00.0095, which has a range of 362 km s"'. 
However, only six sources have total velocity ranges larger than 
100 km s"', and the majority (~60per cent) have total velocity 
ranges of 20kms"' or less. In the lower right panel of Fig.[8]we 
present the isotropic H2O maser luminosity distribution which is 
discussed in the following paragraphs. A summary of minimum, 
maximum, mean, and median values can be found in the last four 
rows of Table[3]for the four distributions presented in Fig.[8] 
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Figure 8. The distributions of derived water maser properties. Tlie up- 
per panels show the peak and integrated flux distribution (the bin size is 
0.5 dex). The lower panels show the total velocity range over which the 
maser spots are seen in each spectrum and the isotropic luminosities. For 
the velocity distribution a bin size of lOkms"' has been used and the pr- 
axis has been truncated at 150 km s" — only one source has a total velocity 
range larger than this. The luminosity distribution uses a bin size of 0.5 dex. 



Kinematic distances have been determined from source veloc- 
ities derived from '^CO JUrquhart et ai]|2008l . |2007|) . or the NH3 
observations p resented here, using the Galactic rotation curve of 
ISrand & Blita ( il993) . assuming the distance to the Galactic cen- 
tre is 8.5 kpc and the angular velocity of the Sun is 220 km s"'. 
We have resolved the twofold kinematic distance ambiguity that 
affects so urces located within the Solar circle using archival HI 
data (see lUrguhart et alJuOl ll for details). These kinematic dis- 
tances have been adjusted for sources found to be associated with 
complexes that have well-determined distances and for sources for 
which maser parallax and photometric distances are available. 

In total we have estimated distances to 573 of the RMS sources 
observed as part of this programme; this includes 295 sources to- 
wards which water maser emission is detected. We use these dis- 
tances to estimate the isotropic luminosity for each water maser by 
integrating all of the emission in e ach channel ab o ve Str and using 
the following equation taken from iAnglada et al.i tl99q) : 



L(H20) 



Lo 



2.30 X 10" 



jSvdV 



Jykms 



D 

kpc 



(7) 



4 RESULTS AND ANALYSIS 

4.1 Detection statistics 

Emission from the ammonia (1,1) inversion transition is detected 
towards 477 sources, with water-maser emission seen towards 
308 sources. These correspond to detection rates of 81 percent 
and 52 per cent (i.e., 477/591 and 308/597) for NH3 and the H2O 
masers, respectively. We find that 273 of the water masers are posi- 
tionally associated with dense gas as traced by the ammonia emis- 
sion. We present a summary of detection rates for the four transi- 
tions for the observed sample in Table[6] as well as a breakdown 
into the four source classifications. The uncertainties in the detec- 
tion rates have been calculated using binomial statistics. 

The detection rates for the Hn/YSO, Hn-region and YSO sam- 
ples are very similar for the ammonia transitions, with the excep- 
tion being NH3 (3,3), which is significantly higher (~60 per cent) 
for the two Hii-region samples compared to the sample consist- 
ing exclusively of YSOs. The NH3 (3,3) inversion transition has a 
higher excitation energy than the NH3 (1,1) and (2,2) lines (~125 K 
above the ground-state energy), and therefore the higher detection 
rate for the NH3 (3,3) transition found towards Hii regions probably 
reflects the presence of warmer gas, which is consistent with these 
regions being more luminous and/or more evolved. 

The water-maser detection rates for the YSOs and Hii-region 
samples are also very similar (~50per cent), indicating that the 
conditions required to pr oduce water-maser e mission are equally 
likely in both phases (cf. lUrquhart et alj|2009r) . The detection rate 
found towards Hii reg ions is comparable with that reported by 
I Kurtz & Hofneii 1 120050 from a set of observations made towards 
a smaller sample conducted with the 100-m Effelsberg telescope 
with a similar sensitivity to the data presented here. 

The ammonia (1,1) and water-maser detection rates are both 
significantly lower (~2 and 3 o", respectively) for the young/old 
classification type, a strong indication that this group includes a 
significant number of evolved stars. The H2O maser and NH3 (1,1) 
and (2,2) detection rates for the young/old group are ~50 per cent 
lower than for the other three source types, while the detection rate 
for the NH3 (3,3) transition is a factor of three lower than the Hii 
sub-samples but only a factor of two lower than the YSOs. 

No ammonia or water maser emission is detected towards ~60 
RMS sources classified as either a YSO, Hii/YSO or Hn region. 
To investigate these non-detections we present a cumulative dis- 
tribution plot showing the heliocentric distance distribution of the 
ammonia and water maser detections and non-detections in Fig.|9] 
Inspection of this plot reveals that the non-detections are pref- 
erentially located at larger distances than sources towards which 
ammonia and/or water masers are detected. For example, we find 
roughly three-quarters of the ammonia and H2O maser detections 
are located within 7 kpc compared with half of the non-detections. 
Given the small beam filling factor estimated in Sect. l3.2l and the 
larger distances associated with the non-detections, we conclude 
that many of the non-detections result from sources falling below 
our observational sensitivity. These non-detections are therefore 
probably due to beam dilution rather than a lack of dense gas asso- 
ciated with any of these sources. 



where D is the distance to the source and the integral extends over 
all components of the spectrum over 3cr. The estimated water maser 
and bolometri c luminosities obtain ed using the bolometric fluxes 
calculated by iMottram et al.l (1201 ih can be found in the last two 
columns of Table|5] The water maser properties derived here will 
be discussed in detail in Sect.l4.5l 



4.2 Detection rates and luminosity 



In Sect. |4.1| we reported the NH3 (1,1) and H2O detection rates for 
the observed sample and as a function of each classification type. 
Overall, we found the NH3 (1,1) and H2O detection rates to be 
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Table 6. Observation and detection statistics. 



NH3 Detection Statistics" 



HtO Detection Statistics 



Source Type Observed Detected Ratio Hyperfine 



Ratio 



(2,2) 



Ratio 



(3,3) 



Ratio 



Observed Detections 



Ratio 



YSO 


275 


235 


0.85±0.02 


199 


Hn region 


214 


167 


0.78±0.03 


129 


Hii/YSO 


71 


60 


0.85±0.04 


48 


Young/old 


31 


17 


0.55±0.09 


10 


Total 


591 


479 


0.81 ±0.02 


387 



0.85±0.03 198 0.84±0.03 

0.77±0.04 149 0.89±0.03 

0.80±0.06 51 0.85±0.05 

0.59±0.16 8 0.47±0.16 



109 


0.46±0.04 


275 


126 


0.75±0.04 


220 


47 


0.78±0.06 


71 


4 


0.24±0.14 


31 



142 


0.52±0.03 


115 


0.52±0.03 


40 


0.56±0.06 


11 


0.35±0.09 



0.81±0.02 406 0.85±0.02 286 0.60±0.02 



597 



308 



0.52±0.02 



" The first, second and third columns give the source type, the number of observations and the total number of sources towards which NH3 emission is detected, 
respectively. In the fourth column we give the overall detection rate. In columns 5-10 we give the total number of sources where the (1,1) hyperfine structure 
is clearly detected and the numbers and detection rates for the (2,2) and (3,3) transitions; the detection ratios for these transitions are given as a proportion 
of the number of NH3 detections (i.e.. Column 3) rather than the number of sources observed. In Columns 1 1-13 we give the number of sources observed for 
H2O masers, the number of detections and the detection rates. 
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Figure 10. Luminosity distribution plots and detection rates as a function of luminosity. In the left and middle panels we present a luminosity histogram and 
cumulative distribution function for the 580 RMS sources with a derived luminosity (black) and the NH3 (1,1) and H2O detections (coloured red and blue 
respectively). In the right panel we present the NH3 (1,1) and H2O detection rates as a function of luminosity; theNHs (1,1) and H2O data points are shown as 
red circles and inverted blue triangles respectively. The error bars presented in these plots have been calculated assuming binomial statistics. The red and blue 
dashed lines show the results of a least-squared fit to the NH3 and H2O data respectively. For the luminosity histogram and cumulative distribution function 
plot we have used a bin size of 0.5 dex, whilst for the detection statistic plot we use a variable bin width so as to include an even number of sources in each 
bin. 




0.0 



5 10 15 

Heliocentric Distance (kpc) 

Figure 9. Cumulative distribution functions showing the relative helio- 
centric distances for the ammonia and water maser detections, and non- 
detections. The distribution of all NH3 and water maser detections is shown 
is shown as a thick black line, while the NH3 (1,1) and H2O and non- 
detections are shown in red, blue and green, respectively. 



~80 per cent and 50 per cent, respectively, and found that these val- 
ues are similar for both the YSO and Hii region samples. In this 
subsection we will use the RMS source luminosities to investigate 



whether or not there is any dependence of the detection rates on 
source luminosity. 

In the left and middle panels of Fig.[TO]we present the lumi- 
nosity distribution for the 580 RMS sources with a derived lumi- 
nosity (thick black line; hereafter we will refer to this sample as the 
RMS luminosity sample) and of the sources towards which NH3 
(red line) and H2O masers (blue line) are detected. In the right panel 
of this figure we present the NH3 (red circles) and water maser (blue 
triangles) detection rates as a function of luminosity. 

Comparing the luminosity distribution RMS luminosity sam- 
ple with that of the NH3 -detected sample reveals them to be indis- 
tinguishable from each other. However, comparing the luminosity 
distribution RMS luminosity sample with that of the H20-detected 
sample we notice a slight increase in the number of more luminous 
sources associated with H2O masers. A Komolgorov-Smirnov (KS) 
test was used to compare the H2O masers sample with the RMS lu- 
minosity sample but found them not to be measurably different. 

The plot presented in the right panel of Fig.llOlreveals the NH3 
(1,1) detection rate to be constant over the whole range of lumi- 
nosities. However, this plot implies that the H2O detection rate has 
a positive dependence on the source luminosity, with the detection 
rate increasing with increasing luminosity. The correlation coeffi- 
cient is ~0.8 which, for A' = 9, is significant with a probability of 
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Figure 11. Comparison plots of the main beam temperatures and FWHM 
line widths of the NH3 (1,1), (2,2) and (3,3) transitions. The dashed red 
lines mark the line of equality. 



<K 1 percent of arising by chance. A linear least-squares fit to the 
data in Fig. 10 predicts the following relationship: 



Log 



RMSh,o 



RMSa 



(0.094 ± 0.014) X log(Lboi) - (0.66 ± 0.07) (8) 



where Lbd is the bolometric luminosity in solar luminosities and 
RMShio/RMSah is the expected detection rate as a function of lu- 
minosity with a sensitivity of ~0.1 Jy beam"' (Icr). The resultant fit 
to the HtO maser detection rate is shown as a dashed blue line in the 
right panel of Fig.[TO] The detections rates of the YSO and Hii re- 
gion samples indicate similar trends, however, lower sample num- 
bers and the larger associated uncertainties prevent us from drawing 
any firm conclusions for these samples. 



4.3 Derived NH3 parameters 

In Table[7]we present mean values for the parameters derived in 
Sect. |3.1| for the whole observed sample and broken down by source 
classification. In this table we include the mean values for the 
'young/old' sample for completeness and, since the values of the 
various parameters associated with the Hii/YSO sample are almost 
indistinguishable from those derived for the Hii region sample, we 
will not discuss the properties of either of these two samples in de- 
tail. We co mplement the resu lts presented in Table|7]with values 
reported bv lPillai et al.l ( l2006r) for a sample of IRDC cores that are 
not coincident with strong compact mid-infrared emission and have 
masses of a few hundred Mq. These IRDC cores therefore consti- 
tute a sample of massive starless cores (as indicated by the lack of 
mid-infrared emission that would indicate the presence of a YSO or 
Hn region) that pre-date the YSO and Hii region phases. Combined, 
these three well-defined classes of objects (IRDC cores, YSOs, and 
Hn regions) should broadly cover the main evolutionary stages of 
massive star formation. 

Looking at the global properties of the whole observed sam- 
ple we find that the mean kinetic temperature is ~20 K, which is 



twice what would be expected from heatin g from the interstel- 
lar radiation field alone (~ 10 K uEvanjI 19991) . These temperatures 
would require the presence of either an internal or external heating 
source, and since we know that the observed objects are associ- 
ated with embedded massive stars, the former is considered more 
likely. We also find the mean line width for the NH3 (1,1) detec- 
tions is ~2kms"', which is approximately ten times broader than 
the thermal linewidth (approximately 0.2 km s"' for gas tempera- 
tures ~20 K). These line widths must be due to systematic motions 
within the gas such as infall, molecular outflows, stellar winds, 
cloud rotation or internal random motions referred to as turbulence. 

Fig.[TT] presents scatter plots comparing the 7,^1, and line 
widths of the NH3 (1,1), (2,2) and (3,3) transitions, revealing a 
strong correlation between the (1,1) and (2,2) line widths and r^b- 
Importantly, the line width ratio for the (1,1) and (2,2) transitions 
is close to unity indicating that the two transitions are tracing a 
similar volume of gas — this was implicitly assumed when deriv- 
ing the physical parameters in Sect. 2.2.1. Comparison of the NH3 
(1,1) and (3,3) transition data reveals a much lower level of correla- 
tion for both parameters, with the NH3 (3,3) transition showing sys- 
tematically lower main beam temperatures and broader line widths. 
The higher gas temperatures required to excite the NH3 (3,3) tran- 
sitions would suggest that the observed emission is emitted from a 
smaller volume of warmer gas located nearer to the embedded ob- 
ject than the NH3 (1,1) emission, which is tracing the cooler outer 
envelope. The broader line widths found for the NH3 (3,3) emission 
are also consistent with this transition tracing warmer and/or more 
turbulent gas closer to the embedded object. 

Inspection of the mean parameters determined for the IRDC 
cores, YSOs, and Hii regions reveals trends for increased linewidth 
and kinetic temperature with evolution phase. Similar trends have 
been reported in the literature and these could be related to the evo- 
lutionary phase of the embedded object as has been suggested by 
ISridharan et al.l ( 120051) from a sample of 56 candidate high-mass 
starless cores, high-mass protostellar objects and UC Hii regions. 
However, as we will show in the following section this interpreta- 
tion is overly simplistic. 

4.4 NH3 correlations 

In the previous subsection we noted a number of trends seen in the 
derived parameters of the IRDC cores, YSO and Hii regions sam- 
ples. We observed a trend for increasing kinetic temperature and 
line width as we moved between the IRDC core, YSO, and Hii re- 
gion samples. Since these three classifications represent an evolu- 
tionary sequence, thes e trends would s eem to suggest an underlying 
connection. However, iMottram et al.l ( 1201 ID showed that there is a 
difference in the luminosity distributions of the YSO and Hii region 
samples, with the Hii regions typically being an order of magnitude 
more luminous. Given the luminosity difference between the sam- 
ples it is possible that the trends seen in the derived parameters are 
actually related to the luminosity of the samples rather being asso- 
ciated with evolutionary progression. 

To investigate the correlation between the derived NH3 pa- 
rameters and source type and/or luminosity, we present six scatter 
plots in Fig.[T2] In the upper panels of this figure we compare the 
sources' column densities with their T^n, AV and Lboi> whilst in 
the lower panels we compare the latter three parameters with each 
other. In these plots we show the distribution of the YSO and Hii re- 
gion sub-samples in red and blue, respectively, and give the Spear- 
man correlation coefficients (rho) in the top left comer of each plot; 
all non-zero correlation coefficients are significant with a probabil- 
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Table 7. Mean values derived from the NH3 transitions. The errors shown in parenthesis are the standard error on the mean. 





Trot 


Tkin 


r 


Log(N(NH3)) 


Tmha.l) 


AV,u, 


rmb(2.2) 


AV(2.2) 


TmbO.i) 


AV(3.3) 


Source Type 


(K) 


(K) 




(cm s-2 


(K) 


(kms-') 


(K) 


(kms-') 


(K) 


(kms-') 


IRDC cores" 


13.34 (0.54) 


13.76(0.56) 


4.25 (0.92) 


15.43(0.16) 


3.12(0.59) 


1.45(0.14) 


1.14(0.41) 


1.51 (0.23) 


1.42(0.39) 




YSO 


17.84 (0.24) 


20.12(0.35) 


2.25(0.11) 


15.22 (0.02) 


1.38 (0.07) 


1.66(0.05) 


0.66 (0.04) 


1.70(0.10) 


0.48 (0.04) 


3.13(0.09) 


Hii region 


20.89 (0.35) 


24.58 (0.57) 


2.33 (0.14) 


15.44 (0.03) 


1.16(0.06) 


2.20(0.11) 


0.64 (0.04) 


2.37(0.17) 


0.50 (0.04) 


4.21 (0.44) 


Hii/YSO 


20.89 (0.43) 


24.41 (0.69) 


1.81 (0.19) 


15.24 (0.05) 


1.11(0.11) 


2.41 (0.21) 


0.64 (0.08) 


2.20 (0.22) 


0.40 (0.07) 


3.28(0.18) 


Young/old 


15.74 (0.59) 


17.18(0.75) 


2.91 (0.61) 


15.29(0.12) 


1.08 (0.22) 


0.74(0.11) 


0.52(0.14) 


1.41 (0.19) 


0.31 (0.13) 


2.47 (0.47) 




19.22 (0.20) 


22.12(0.30) 


2.24 (0.08) 


15.30 (0.02) 


1.27(0.04) 


1.91 (0.06) 


0.65 (0.03) 


2.00 (0.08) 


0.47 (0.03) 


3.62 (0.20) 


" The values presented in this 


row have been taken fromlPillai et alj (20061) 
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Figure 12. Scatter plots showing the correlation between the column density, Tkin, AV and Lboi. The YSO and Hn region sub-samples are shown in red and blue 
respectively. The Spearman correlation coefficients are given in the top left corner of each plot. The black dotted line shows the result of linear least-squared fits 
to all of the plotted data, whilst the red and blue dashed lines show the results of linear least-squared fits to die YSO and Hn region sub-samples respectively. 



ity of arising by chance of <k 1 per cent. We find a correlation be- 
tween five of the six pairs of plotted parameters and in these cases 
we have performed a finear fit to the combined YSO and Hit re- 
gion samples which is shown as dotted black line. In cases where 
we find a moderate correlation (i.e., 0.5<rho <0.75) we have made 
separate linear fits to the YSO and Hit region sub-samples to inves- 
tigate possible differences in the overall trends; these are shown as 
dashed red and blue lines respectively. The intercepts and gradients 
for all of these fits are presented in Table[8] 

Visual inspection of the upper panels of Fig.[T2]reveal there is 
no correlation between column density and Ty,, and only a weak 
correlation between column density and the other two parameters. 
The correlation between the Tun, AV and Lboi parameters is signifi- 
cantly stronger than found for the column density with values of rho 
> 0.5. However, we note that the correlation for these parameters is 
lower for the YSO and Hit region samples with only a weak corre- 



lation being found for the Hn regions. The correlation between the 
Tkin and AV (lower left panel of Fig.ll2t is particularly strong with 
a correlation coefficient of 0.66. Moreover, the YSOs and Hit re- 
gions dominate either end of what appears to be fairly continuous 
distribution. This would seem to support the trend of increasing 
line width and kinetic temperature with evolution touched upon in 
Section |473l However, turning our attention to the lower middle and 
right panels of Fig.[T2]we find that both of these parameters are also 
correlated with the sources' bolometric luminosity. 



IChurchwell et al.l ( 1 19901) found a similar, although weaker, 
correlation between the kinetic temperature and linewidth with far- 
infrared luminosity from NH3 observations of a smaller sample of 
UC Hit regions (84) made with the Effelsberg 100-m telescope. The 
authors suggested that the line width-luminosity correlation can be 
understood as either: 1) the amount of mechanical energy deposited 
into the ambient molecular cloud (e.g., via molecular outflows, stel- 
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Table 8. Gradients and intercepts of linear least-squared best fits to the plots 
presented in Fig. ll2l These fits have been performed in log-log space. 



Plot 


Group 


Intercept 


Gradient 


rho 


N(NH3) vs. AV 


All 


14.28±0.31 


4.10±1.16 


0.46 


N(NH3) vs. Lboi 


All 


11.62±().25 


0.93±0.06 


0.37 




All 


0.71±0.03 


0.15±0.01 


0.51 


Ttin VS Lbol 


YSO 


0.69±0.04 


0.16±0.01 


0.49 




Hii 


0.21±0.15 


0.26±0.03 


0.30 




All 


10.59±0.07 


4.21 ±0.05 


0.66 


rk,„ vs. AV 


YSO 


10.69±0.07 


4.04±0.05 


0.63 




Hi: 


13.55±0.32 


3.29±0.15 


0.49 




All 


-1.10±0.06 


0.34±0.02 


0.63 


AV vs. Lbo, 


YSO 


-1.28±0.11 


0.40±0.03 


0.56 




Hi: 


-1.71±0.25 


0.46±0.06 


0.44 



lar winds) being directly proportional to the luminosity of the em- 
bedded source; or 2) the line width-clump mass relationship (i.e., 
iLarsonlligSlh . with the more luminous, and presumably more mas- 
sive, stars forming out of more massive clumps. The correlation be- 
tween kinetic temperature and luminosity would indicate that the 
heating of the gas is dominated by the radiation deposited by the 
embedded star into its natal environment. 

The distributions presented in the right panels of Fig.[T2] sug- 
gest that there is a dependence of line width and kinetic temperature 
on bolometric luminosity. It is clear from these correlation plots 
that the embedded sources are having a measurable effect on their 
surrounding molecular environment. Given that these observations 
are sampling clumps with sizes ranging from 0.1-0.5 pc it is reason- 
able to conclude that it is the energy output of a young massive star 
that is driving the increases in temperature and line width. However, 
it is not yet clear if there is also a weaker underlying dependence 
on the evolutionary phase of the embedded sources. 

Inspection of the plots and fit parameters presented in Table[8] 
would suggest there is a significant difference between the gas tem- 
perature and line width, and the gas temperature and the luminosity 
for the Hii region and YSO sub-samples. There is no significant dif- 
ference between the line width-luminosity relationship for the two 
samples, which would suggest that it is differences in the kinematic 
gas temperature that is the parameter of interest. In an attempt to 
investigate the gas temperature and line width differences for the 
Hn regions and YSO sub-samples we have removed the luminos- 
ity dependence using their respective fit parameters; this effectively 
normalises both parameters. In Fig.[T3]we present normalised his- 
tograms and cumulative distribution plots showing the differences 
between the YSO and Hn region samples after the luminosity de- 
pendence has been removed. These plots reveal little difference be- 
tween the two normalised parameters. Komolgorov-Smirnov (KS) 
tests used to compare the two samples found that the probabili- 
ties that YSOs and Hn regions were drawn from the same parent 
population were 3 per cent and 4 per cent for the normalised kinetic 
temperature and line width parameters respectively. We therefore 
conclude that the observed trends in Table 7 are more likely to be 
due to the energy output of the central source and/or the line width- 
clump mass relationship than the evolutionary state of the embed- 
ded source. 



4.5 H2O maser correlations 

Given the ubiquitous nature of water masers within the Galactic 
plane it is possible that a number of the detected masers are the re- 
sult of chance line of sight alignments within the field or a bright 
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Figure 13. Histograms and cumulative distribution plots showing the dif- 
ferences between the YSO and Hii region samples after the luminosity de- 
pendence has been removed. The distributions of the YSO and Hii regions 
samples are shown in red and blue respectively. 
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Figure 14. Distribution of the differences between the velocity of the bright- 
est water maser spot and that of the ammonia emission. 



maser located within a telescope sidelobe rather than genuine asso- 
ciation; however, given the GBTs off-axis design, sidelobe contam- 
ination is significantly lower than for traditional single-dish tele- 
scopes and so the latter possibility is less of a concen£|. In this 
section we will investigate the degree of correlation between the 
water maser and NH3 emission in order to ascertain whether they 
are arising from the same molecular gas, and if so, whether there is 
a direct connection between the maser and the embedded massive 
young star. 

In Fig.[T4] we present a histogram of the difference in veloc- 
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ity between the most intense water maser spot and that of the dense 
molecular gas traced by the NH3 emission. This plot reveals that the 
velocities of the two transitions are in excellent agreement, with a 
Pearson correlation coefficient of 0.92. The mean velocity diifer- 
ence is -3.8 with a standard deviation of ~20kms"'. Computing 
the standard error on the mean, we find that this small, negative 
offset from zero is significant at the > 3-cr level. However, the stan- 
dard deviation is inflated by a small number of high-velocity out- 
liers in the distribution. If we exclude the 25 sources whose the ve- 
locity difference is larger than 20 km s"' , we find that the mean and 
standard deviation reduce to ~ -1.1 kms"' and 6kms"', respec- 
tively. The associated standard error on the mean is ~ 0.38 km s"', 
so the revised mean is still nearly 3 a away from zero. Thus there 
appears to be a real offset from zero in the mean relative veloc- 
ity of the maser emission. The distribution also has a slight skew 
(~0.26kms"') towards blue-shifted sources, with a negative me- 
dian value (-0.82 km s'). So not only are there more blue-shifted 
than red-shifted masers but they also have higher relative veloci- 
ties. This probably means that the blue-shifted maser spots tend to 
be slightly brighter than their red-shifted counterparts and so more 
of them are detected. 

In Fig.[T5]we plot the NH3 velocities against the velocity of the 
brightest maser spot. The black horizontal lines shown in this plot 
indicate the velocity range over which maser spots are seen towards 
each source. The data presented in Fig.[T5]illustrate the correlation 
between the molecular gas and the range of maser velocities and 
the excellent correlation between the velocity of the brightest maser 
spots and the molecular gas. For the vast majority of sources, the 
peak maser velocity and NH3 velocity are within 20 km s"' of each 
other, but there are 25 sources in which the velocity difference is 
larger. The light yellow stripe running through the plot presented in 
Fig.[T5]indicates the region where the H2O maser velocity is within 
±20kms"' of the NH3 emission. This number is further reduced 
if we take the total velocity range of the maser spots seen towards 
each source into account; of these 25 sources we find that the NH3 
velocity of five lie within the maser velocity range, and five others 
are located within 20 km s"' of a maser spot. 

There are 15 sources (~6 per cent of our sample) for which 
there is >20kms"' difference between the velocity of the nearest 
maser spot and the velocity of the m olecular cloud. Ta r geted NH3 
and H2O observation s conducted by IChurchwell et al.l i ll 9901) and 
lAnglada et al.l i ll 9961) (that have approximately the same resolution 
and sensitivity to the observations we present here) reported the de- 
tection of a similar propo rtion of sources with large velocity offsets. 
IChurchwell et al.l jl990l) concluded that these are simply serendipi- 
tously matched sources that are located nearby on the sky but which 
arise from physically distinct regions located along the same line of 
sight. However, the strong variability of water masers and/or lim- 
ited sensitivity could also explain the large velocity differences. 

In a recent study of water masers by iBreen et alj | |2010|) . 
large changes in the velocity of the brightest maser peak were 
found for four sources by comparing the emission spectra taken 
in two epochs separated by a period of 10 months. Taking source 
G336. 983-0. 183 as an example, they found only a single maser 
spot at approximately -75 km s"' in their 2003 observations. How- 
ever, in the spectrum obtained in 2004 they found two maser spots 
at approximately -75 and 45 km s"'. It is therefore possible that the 
strong variability of water masers can sometimes result in the large 
velocity differences observed towards some sources. So although 
we find large velocity differences between a small number of wa- 
ter masers and the molecular gas it is not yet clear that these arise 
from physically separate regions. It is clear, however, that there is 
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Figure 15. Compaiison of the velocities of the NH3 emission and H2O 
maser. The red circles show the velocity of the NH3 emission and the ve- 
locity of the brightest H2O maser spot, whilst the horizontal black line in- 
dicates the velocity range over which the maser emission is detected. The 
dashed line shows the line of equality where the two velocities are equal and 
the the light yellow stripe mnning through the plot indicates region where 
the H2O maser velocity is within ±20kms"' of the NH3 emission. 



a strong correlation between velocities of the NH3 and H2O tran- 
sitions for ~95 percent of the sources and it is therefore likely that 
they are associated with each other. 

By making the assumption that all of the detected water 
masers are associated with the molecular clouds, we can inves- 
tigate the distribution of high-velocity features, which we define 
as those which vary from the clouds' syste mic velocity by more 
than 30km S-' (c.f ICaswell & Breenll20lol) . We find 10 sources 
that show both high-velocity red- and blue-shifted features, 33 that 
show only blue-shifted features, 17 that show only red-shifted fea- 
tures and 213 sources that show no high-velocity features. We note 
that there are almost twice as many sources with high-velocity blue- 
shifted features. A similar distribution o f high-velocity red- and 
blue-shifted features has been reported bv iCaswell & BreenI l|2010|) 
who interpreted the excess of blue-shifted features as a sign of the 
youth of the embedded object as there was also a strong correla- 
tion with these blue-shifted features and methanol masers, which 
are considered to be exclusively associated with the earliest phases 
of high-mass star formation. 

This excess in blue-shifted, high-velocity maser emission re- 
flects the small but significant skew in the general distribution de- 
scribed above and is likely to be related to expansion or contrac- 
tion in the core or envelope gas. Further analysis would require 
detailed modelling of the velocity field and density distribution. 
We compared the maser parameters with the column densities and 
line widths obtained from the ammonia data and as a function of 
source type to identify any correlations that might provide any fur- 
ther insight into this excess of blue-shifted emission, however, these 
proved inconclusive. 

Having established that the observed water-maser emission is 
very likely to be associated with the same molecular clouds within 
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Figure 16. Comparison of the bolometric luminosity of RMS sources with 
the isotropic water maser luminosity. The blue dots show the luminosities 
of all maser detections for which we have resolved the distance ambiguity 
while the red dots show the source luminosities of a distance limited sample 
(see text for details). The dashed line shows the linear bisector least-squares 
fit to the data. The partial Spearman correlation coefficient and the and re- 
sults of the fit to the data are given in the top left of the plot. 



which the MYSOs and Hii regions are forming, the next step is to 
see if we can identify a direct connection between the masers and 
the embedded young massive star. In Fig.[T6]we present a scatter 
plot comparing the isotropic maser luminosity calculated in Sec- 
tion |3.4| with the RMS luminosity. On this plot we show the distri- 
bution of all of the H2O detections in blue and a smaller distance 
limited sample is shown in red. The distance limited sample in- 
cludes all sources with heliocentric distances less than 3 kpc with 
luminosities greater than 1000 Lq and should therefore be relatively 
unaffected by the Malmquist bias that can lead to false correlations. 
The dashed lines show the bisector least-squared fits to the data and 
clearly shows that the water maser luminosity is correlated with the 
bolometric luminosity for the whole and distance limited samples, 
and that the fits are in good agreement with each other. As an ad- 
ditional check we calculate the partial Spearman correlation coeffi- 
cient using the square of the distance as an independent parameter; 
a value of 0.44 would suggest the correlation is real. From the lin- 
ear bisector least-squares fit to all of the plotted data we obtain the 
following relationship: 



LH20 - a(L\,oi) 



(9) 



where Lboi is the RMS bolometric luminosity and a = (7.1 ± 0.3 x 
10"'^ and b = (1.47 + 0.76). We conclude from the correlation of 
the cloud and water-maser velocities and the bolometric and maser 
luminosity that there is a strong dynamical relationship between 
the embedded young massive star and the H2O maser as has been 
previ ously suggested from s i milar studies of star forming regions 
(e.g.. lChurchwell et al.lll990l : lAnglada et aljl996h . 



4.6 Ammonia (3,3) maser? 

We detected thermal NH3 (3,3) emission towards 290 sources 
which corresponds to approximately 50 per cent of the sample ob- 
served. In all but one of these cases we found that the observed 
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Figure 17. The NH3 (1,1), (2,2) and (3,3) spectra detected towards 
G030.7206-00.0826. The NH3 (3,3) emission is clearly made up from two 
components, a broad component thermal in nature and a stronger thin com- 
ponent possibly indicating the presence of an ammonia maser; the fits to 
these two components are coloured red and blue respectively. 



emission could be approximated by a single Gaussian profile. How- 
ever, the NH3 (3,3) emission detected towards the MSX source 
G030.7206-00.0826 consists of two distinct components, a broad 
component associated with the thermal emission from the embed- 
ded source and a second narrow-lined component. The ammonia 
spectra for this source are presented in Fig.fTTl 

We have fitted the NH3 (3,3) emission with two Gaussian pro- 
files to account for the broad and narrow-lined component. The 
fitted linewidth of the narrow component is 1.1 kms"' and it has 
a peak flux density of ~5 Jy, which makes it the narrowest NH3 
(3,3) line detected in the sample. Since this narrow-lined compo- 
nent does not have a corresponding feature seen at the same ve- 
locity in the NH3 (1,1) and (2,2) emission profiles it is unlikely 
to be thermal in nature. Given its unusually narrow line width this 
is likely to be a maser transition. It is not possible to confirm the 
nature of this emission feature with the current data and we await 
further observations. 

Another susp ected NH3 maser ( G23. 33-0.30) has recently 
been reported by IWalsh et al.l ( 1201 ih from the H2O Souther n 
Galactic Plane Survey (HOPS) 12-mm survey jWalsh et al.ll2008l) . 
G23. 33-0.30 has a flux density of 9.7 Jy, which is comparable with 
the flux den s ity of the narrow-lined emission feature reported here. 
IWalsh et all ( l201lh note that although the NH3 (3,3) transition has 
been detected as a maser in a number of star forming regions (e.g., 
PR 2 1 (OH) — iMaiigum & WoottenI 1 1 9941 ; W51 — IZhangetal] 
I1995I ; G5.89-0.39 — iHunter et al.l2008h it has never been stronger 
than 0.5 Jy and thus, if these two masers can be verified they would 
be the strongest NH3 (3,3) masers reported to date. 

G030.7206-00.0826 is associated with a b right unresolved ra- 
dio source (~0.5 Jy at 5 GHz; IWhite et alj2005h and has been clas- 
sified as an Hn regi on with a bolometric luminosity of ~ 4 x lO'* Lq. 
IWalsh et al.l ( I2OIII) found G23. 33-0.30 to be positionally coinci- 
dent with an IRDC and a cold dust continuum source and thus 
it therefore seems likely that these maser candidates arise in star 
forming regions. However, since only two NH3 (3,3) maser can- 
didates have been reported from the 100 square degrees surveyed 
by HOPS and from the GBT observations of ~600 high-mass 
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star forming regions reported here, we would conclude that bright 
masers of this type are relatively rare. 



5 SUMMARY AND CONCLUSIONS 

We have used the 100-m Green Bank Telescope to conduct a pro- 
gramme of 22-24 GHz spectral line observations towards 597 mas- 
sive young stellar objects (MYSOs) and compact and ultra-compact 
Hn regions identified by the Red MSX Source (RMS) survey. We 
targeted the 22 GHz water (H2O) maser transition and thermal am- 
monia (NH3) (1,1), (2,2) and (3,3) inversion transitions. These ob- 
servations have an angular resolution of ~30" and typical Act sen- 
sitivity of r,^i, = 0.2 K per 0.32-km s"' channel. 

We detect water maser emission towards 308 RMS sources 
with an overall detection rate of ~50 per cent. NH3 emission is de- 
tected towards 479 sources (~80 per cent of the sample) with the 
NH3 (1,1) hyperfine structure being clearly detected towards ~400 
sources, and thus allowing the kinetic temperature, optical depth 
and column density to be determined for a large number of massive 
star-forming regions. The average kinetic temperature, FWHM line 
width and column density for the sample are approximately 22 K, 
2kms"' and 2 x 10'^ cm"-, respectively. 

We find the detection rates for both the HtO masers and NH3 
emission is similar for both MYSO and Hn region sub-samples. 
Combining these results with a large database of complementary 
observations we analyse these parameters as a function of luminos- 
ity and evolutionary phase, and correlate these physical conditions 
with HtO maser emission. Our main findings are as follows: 

(i) There is no significant difference in the H2O maser detec- 
tion rate for Hii regions and MYSOs which would suggest that the 
conditions required to produce maser emission are equally likely 
in both phases. Comparing the detection rates as a function of lu- 
minosity we find the H2O detection rate has a positive dependence 
on the source luminosity, with the detection rate increasing with 
increasing luminosity. 

(ii) We find that the NH3 (1,1) line width and kinetic temper- 
ature are correlated with luminosity. It is clear from these corre- 
lations that the embedded sources are having measurable effects 
on their surrounding molecular environments. Finding no underly- 
ing dependence of these parameters on the evolutionary phase of 
the embedded sources, we conclude that the observed trends in the 
derived parameters are related to the energy output of the central 
source and/or the linewidth-clump mass relationship. 

(iii) The velocities of the peak H2O masers and the NH3 emis- 
sion for each relevant source are in excellent agreement with each 
other with a correlation coefficient of 0.92, which would strongly 
suggest an association between the dense gas and the maser emis- 
sion. Moreover, we find the bolometric luminosity of the embedded 
source and the isotropic luminosity of the H2O maser are also cor- 
related. We conclude that there is a strong dynamical relationship 
between the embedded young massive star and the H2O maser. 

(iv) We find excellent agreement between velocities of sources 
using ''CO observations and those obtained from NH3 transitions. 
We find only three sources where the velocity assigned using the 
CO data is incorrect. 

These observations are a first step in examining the global 
characteristics of this Galaxy-wide sample of massive young stars 
and will form the cornerstone for more detailed studies of well se- 
lected sub-samples. 
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APPENDIX A: COMPLEX AMMONIA SPECTRA 
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Figure Al. Complex spectra obtain towards six RMS sources that show evidence of both emission and absorption. 



